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Abstract 

In stars with temperatures above 20- 10 6 K, hydrogen burning is dominated by the 
CNO cycle. Its rate is determined by the slowest process, the 14 N(p,7) 15 reaction. 
Deep underground in Italy's Gran Sasso laboratory, at the LUNA 400 kV acceler- 
ator, the cross section of this reaction has been measured at energies much lower 
than ever achieved before. Using a windowless gas target and a 4ir BGO summing 
detector, direct cross section data has been obtained down to 70keV, reaching a 
value of 0.24 picobarn. The Gamow peak has been covered by experimental data for 
several scenarios of stable and explosive hydrogen burning. In addition, the strength 
of the 259 keV resonance has been remeasured. The thermonuclear reaction rate has 
been calculated for temperatures 90-300 -10 6 K, for the first time with negligible 
impact from extrapolations. 
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1 Introduction 



Charged particle induced thermonuclear reactions in a star take place in a 
narrow energy window called the Gamow peak, far below the Coulomb barrier. 
Their cross section cr(E) drops steeply with decreasing energy and can be 
parameterized as [1]: 

a(E) = ^le~^ (1) 

where S(E) is the astrophysical S-factor, and rj oc E~ - 5 is the Sommerfeld 
parameter 1 . The very low value of cr(E) at stellar Gamow peak energies pre- 
vents a measurement in a laboratory at the surface of the Earth, where the 
signal to background ratio is too small because of cosmic ray interactions. 
In order to overcome this difficulty, a 50 kV [2] and a 400 kV [3] accelerator 
have been installed in the Laboratory for Underground Nuclear Astrophysics 
(LUNA), deep underground in Italy's Gran Sasso laboratory. At LUNA, two 
reactions from the hydrogen burning p-p chain have been studied for the first 
time directly at their respective solar Gamow peak: 3 He( 3 He,2p) 4 He [4,5] and 
2 H(p,7) 3 He [6]. 

The second process of stable hydrogen burning in stars, the carbon-nitrogen- 
oxygen (CNO) cycle, dominates for 20<T 6 < 130 (T 6 : stellar temperature in 
10 6 K), corresponding to Gamow peak energies of 30-110keV. The slowest 
process, the 14 N(p,7) 15 reaction, has been studied previously [7, and ref- 
erences therein], but only one work reported data at these low energies [8]. 



+39 049 827 7145. 

1 In the present work, E denotes the energy in the center of mass system, and E p 
is the projectile energy in the laboratory system. 
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Reaction rate compilations [9,10,11] have been based mainly on ref. [7], with 
data down to E— 181 keV, and on ref. [8]. Subsequently, the extrapolation to 
solar energy of ref. [7] for capture to the ground state in 15 has been chal- 
lenged by several works, on theoretical [12] and indirect [13,14,15] grounds. 
These considerations [12,13,14,15] were then experimentally confirmed at E — 
119-367keV [16,17], indicating an extrapolated S-factor at solar energy that 
is a factor 2 smaller than the value adopted in the compilations. Recently, 
an independent study reported cross section data for E— 134-482keV, with 
identical conclusions for the total extrapolated S-factor [18]. 

In the present work, the results of a novel 14 N(p,7) 15 experiment are pre- 
sented. The previous studies (table 1) relied on extrapolations in order to 
predict the S-factor at astrophysical energies [7,12,13,14,15,16,17,18] or cover 
only a narrow energy range [8]. The present work follows a different approach. 
Cross sections are measured directly, with high statistics and over a wide, 
astrophysically relevant energy range. 

2 Experiment 

Equation 1 predicts a factor 200 drop in the yield from E— 119 keV [16,17] to 
70keV (the aim of the present study), so a new setup with a thin, pure gas 
target and a high-efficiency annular BGO detector had to be developed [19]. 
In the experiment, the proton beam with E p = 80 - 250 keV and 0.3 mA typical 
current from the LUNA 400 kV accelerator [3] passed several apertures (final 
aperture: 40 mm long, 7mm diameter) and entered a windowless gas target. 
The 12 cm long target cell is placed inside the borehole (6 cm diameter) of a 
28 cm long annular bismuth germanate (BGO) crystal in 4n geometry, ensuring 
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^70% peak detection efficiency for 7MeV 7 rays [19]. The beam is stopped 
on a copper disk that also serves as the hot side of a calorimeter with constant 
temperature gradient [19]. 

The beam energy has been calibrated previously for E p = 130-400keV [3], 
and the calibration is extrapolated down to £ , p = 80keV, resulting in 0.3 keV 
uncertainty. The target gas was 1 mbar nitrogen (corresponding to typically 
lOkeV target thickness), 99.9995% chemically pure. The pressure was mon- 
itored with a capacitance gauge and kept constant to better than 0.25 % by 
a feedback system. The beam heating effect [20] reduced the target density 
by up to 15% with incident ion beam. Using the resonance scan technique 
[1,20], this correction was measured for the present setup, prior to the actual 
experiment. 

Because of the high absolute detection efficiency and of the near 4ir geometry 
of the BGO detector, 7 rays emitted in a cascade are with high probability 
summed into a peak at E 1 = Q + E (Q — 7.297 MeV is the Q value for the 
14 N(p,7) 15 reaction) in the spectrum (fig. 1). Next to the summing peak are 
unresolved lines at 6.172 and 6.792 MeV, the energies of two secondary 7 rays. 
A peak at 5.6 MeV results from the 2 H(p,7) 3 He beam induced background 
reaction [21]. The broad structure at 12 MeV is caused by the 15 N(p,7) 16 re- 
action (the target gas has natural isotopic composition, 0.4 % 15 N), with a con- 
tribution from the 11 B(p,7) 12 C reaction. The latter reaction also gives 7 rays at 
16 MeV. The beam induced background in the region of interest (ROI) results 
from the Compton continuum due to high energy 7 rays (from 15 N(p,7) 16 
and n B(p,7) 12 C), and from the 13 C(p,7) 14 N reaction (Q = 7.551 MeV). The 
number of events from 13 C(p,7) 14 N (E 1 xs 7.7 MeV) has been evaluated by tak- 
ing monitor 7 spectra with helium instead of nitrogen as target gas. The bulk 
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of the 13 C (due to hydrocarbons from pump oil) seen by the beam has been 
localized at the beam stop [21]. The n B has been localized on the final collima- 
tor [21], and it is due to impurities in the collimator material. The laboratory 
background in the ROI is due to (11,7) reactions (caused by (a,n) neutrons 
from the laboratory background) mainly in the stainless steel target chamber, 
and it has a constant counting rate [21]. The total 7 ray background amounts 
to less than 20 % of the counts observed in the ROI for each run, except for 
the two lowest energy points: At E = 80 keV, the background is 45 % of the 
collected events. At E=70keV, with 53 days of running time, there are 11 
counts/day from the reaction, 22 counts/day from laboratory background and 
1 count/day from beam induced background. 

The 7 ray detection efficiency depends on the measured, energy dependent 
branching ratios [17], and for £<119keV, where there is no data, on R- 
matrix extrapolations [17]. However, this dependence is weakened [22] by the 
flatness of the efficiency curve [19] and the dominance of the 6.792 Me V tran- 
sition [7,12,14,16,17,18,23]. As a result, the branching ratios contribute 3% 
systematic uncertainty in the detection efficiency, with an additional 1 % from 
the modeling of the detector [19] and 1.5% from the absolute efficiency cal- 
ibration, giving a total of 3.5%. The angular distributions for primary and 
secondary 7 rays from the five strongest transitions have been measured pre- 
viously above and below the 259 keV resonance [17], and isotropy was found 
for all 7 rays except for the primary from capture into the state at 6.792 MeV. 
Since both the 6.792 MeV secondary 7 ray and the sum peak are fully included 
in the ROI, the detection efficiency is unaffected [22]. Details on experiment 
and analysis can be found elsewhere [24,25] and will be published separately. 
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3 Results 



The present S-factor data (fig. 2) reach a much lower energy than the previous 
direct experiments [7,8,16,17,18], while overlapping over a wide energy range. 
The statistical uncertainties are lower, and the systematic uncertainties are 
comparable or lower (table 1). The systematic uncertainty for the lowest en- 
ergy point in the present study is 7%, dominated by beam energy calibration 
(5%), detection efficiency (3.5%) and effective target density (3%). 

For the lowest energy points of refs. [16,17], not all transitions given in table 
1 were measured, and R-matrix extrapolations [17] had to be added, giving a 
7 % contribution. An analogous procedure had to be applied to the lowest data 
points of ref. [18], with R-matrix fit results [18,17] contributing 26% to the 
lowest data point. In the resulting total S-factor picture (fig. 2), the present 
data are systematically lower than ref. [8] , whereas good agreement is obtained 
with ref. [7] in the overlapping energy region. Excellent agreement is obtained 
between the present data and refs. [16,17]. Good agreement is observed with 
ref. [18], except for E = 185-215keV, where ref. [18] is systematically lower 
than both the present data and refs. [16,17]. 

In addition to the cross section measurement, the strength w-f of the 259 keV 
resonance has been determined: tvy = 12.8±0.3 sta t±0.4 syst meV, in excellent 
agreement with previous works: 12.9±0.4 stat ±0.8 syst meV [17], 13.5±1.2meV 
[18] and 14±1 meV [26]. The value by ref. [26] has been adopted by the Nuclear 
Astrophysics Compilation of Reaction Rates (NACRE) [11]. 

The present astrophysical S-factors, corrected for the electron screening effect 
[27] in the adiabatic limit (10 % correction for E = 70 keV, 3 % for E = 150 keV), 
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and the present uj~/ value have been used to compute the thermonuclear re- 
action rate (fig. 3). For energies i?<70keV, the S-factor has been assumed 
to be constant and equal to the 70keV value. For temperatures T 6 >90, the 
experimental data from the present work account for more than 90 % of the 
area under the Gamow peak, with the remaining 10% depending on the as- 
sumption made for E < 70 keV. For 90 > T & > 60, data account for 50 - 90 % 
of the area under the Gamow peak. 

For T e > 180, the rate is dominated by the 259 keV resonance, and the present, 
lower resonance strength leads to a rate that is systematically 10% lower 
than NACRE. For T 6 < 180, nonresonant capture becomes more and more 
important, and the present, experiment-based rate is up to 40 % lower than 
NACRE. 

4 Astrophysical consequences 

In several astrophysical scenarios the present, revised reaction rate for T 6 > 60 
has direct consequences. After the end of their helium burning phase, low- 
mass stars burn hydrogen and helium, respectively, in two shells surrounding 
a degenerate carbon-oxygen core and give rise to the so-called asymptotic gi- 
ant branch (AGB) [28] in the Hertzsprung-Russell diagram. Flashes of helium 
burning spawn convective mixing in a process called dredge-up, transporting 
the products of nuclear burning from inner regions of the star to its surface, 
where they are in principle accessible to astronomical observations. The tem- 
perature in the hydrogen burning shell is of the order of T 6 = 50 - 80 for a 2 M & 
AGB star (M : mass of the Sun) with metallicity Z = 0.01 [29]. It has been 
shown that an arbitrary 25% reduction of the 14 N(p,7) 15 rate with respect 
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to NACRE leads to twice as efficient dredge-up of carbon to the surface of the 
star [29] . At these temperatures, the present rate is 40 % below NACRE (fig. 
3). Recently, a simulation for a 5M , Z = 0.02 AGB star [30] found stronger 
thermal flashes for a reduced CNO rate, consistent with the result by ref. [29]. 
Explosive burning in novae [31] takes place at even higher temperatures, typ- 
ically T 6 ^200, through the hot CNO cycle. The abundance of 15 N (daughter 
of 15 0) in nova ashes depends sensitively on the 14 N(p,7) 15 rate [32]; the 
present, significantly more precise rate reduces the nuclear uncertainty of the 
predicted abundance. 



5 Summary 



The total cross section of the 14 N(p,7) 15 reaction, the bottleneck of the 
CNO cycle, has been measured for E — 70- 228 keV, with typically 3 % (at 
most 10 %) statistical and 5 % (at most 7 %) systematic uncertainty. For the 
first time, precision cross section data has been obtained directly at energies of 
hydrogen burning in AGB stars. The strength of the E R = 259 keV resonance 
has been determined with improved precision. The thermonuclear reaction 
rate for several scenarios of stable and explosive hydrogen burning has been 
calculated directly from the present experimental data, with negligible impact 
from assumptions made for lower energies. Several significant consequences 
of the present, experiment-based rate for the evolution of AGB stars and for 
nucleosynthesis in novae have been discussed. 
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Table 1 

Cross section measurements of 14 N(p,7) 15 at low energy. Capture to the states 
at 6.859 and 7.276 MeV in 15 is negligible [7,14]. Uncertainties given refer to the 
lowest energy point. T 6 min is the lowest temperature where more than 90 % of the 
integral over the Gamow peak is covered by experimental data. 



Year 


Ref. 


Data on capture to states: 


E [keV] 


Uncertainty 
stat. syst. 


min 
J 6 


1957 


[8] 


All (activity measurement) 


93 - 126 


46 % 15 % 


none 2 


1987 


[7] 


GS 5.181 5.241 6.172 6.792 


181 - 3600 


32 % 13 % 


180 


2004/5 


[16,17] 


GS 5.181 5.241 6.172 6.792 


119 - 367 


7% 10% 


150 


2005 


[18] 


GS 6.172 6.792 


134 - 482 


38 % 9.4 % 


170 


Present work 


All (47r summing crystal) 


70 - 228, 259 


10% 7% 


90 



2 The energy range of ref. [8] is so narrow that it does not cover more than 50 % of 
the Gamow peak for any given temperature. 
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Fig. 1. Typical 7 ray spectrum in the 4-7T BGO detector, at E = 90 keV, livetime 11.4 
days, accumulated charge 231 C. The counts from the reaction to be studied are 
shaded in dark gray. The laboratory background [21], normalized to equal livetime, 
is shaded in light gray. The most important components of background induced by 
the ion beam are indicated. Inset: Region of interest (ROI) for the present study, 
with the beam induced background in the ROI indicated by the black filled area 
(see text). 
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Fig. 2. Astrophysical S-factor for the 14 N(p,7) 15 reaction from the present work 
(filled squares) and from previous studies: circles [8], inverted triangles [7], diamonds 
[16,17], triangles [18]. Error bars (±lcr statistical uncertainty) are only shown where 
they are larger than the symbols used. The Gamow peak for Tq = 80 is also shown. 
The systematic uncertainties are given in the text and in table 1. 
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Fig. 3. Thermonuclear reaction rate relative to the NACRE [11] rate. Dot-dashed 
line: NACRE rate. Solid (dashed) line: present work and more than 90 % (more than 
50 %) of the Gamow peak covered by experimental data. Dotted line: Extrapola- 
tion-based rate from ref. [17]. The shaded areas indicate quoted upper and lower 
limit for the NACRE rate [11] and ±1<t statistical uncertainty for the rate from ref. 
[17] and the present work. 
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